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Nuclear /3-decay in magic nuclei is investigated, taking into account the coupling between particle 
and collective vibrations, on top of self-consistent random phase approximation calculations based 
on Skyrme density functionals. The low-lying Gamow-Teller strength is shifted downwards and at 
times becomes fragmented; as a consequence, the /3-decay half-lives are reduced due to the increase 
of the phase space available for the decay. In some cases, this leads to a very good agreement between 
theoretical and experimental lifetimes: this happens, in particular, in the case of the Skyrme force 
SkM* *, that can also reproduce the line shape of the high energy Gamow-Teller resonance as it was 
previously shown. 

PACS numbers: 21.60.Jz, 23.40.-s, 26.30.-k 


Weak interaction processes involving atomic nuclei 
have been the object of continuous interest for many 
decades [I]. The simplest process, that is, nuclear ,8- 
decay, is mainly determined by the allowed Gamow-Teller 
(GT) type of transition M. As a rule, most of the GT 
strength is concentrated in the high-energy GT resonance 
(GTR); consequently, the fraction that contributes to the 
8-decay, being in the energetically allowed region called 
“8-decay window”, is small. The distribution of the GT 
strength is governed both by the nuclear shell structure 
(in particular, by the spin-orbit splitting between neu¬ 
tron and proton states), and by the still underconstrained 
spin-isospin channel of the nuclear effective interaction, 
so that 8-decay and GTR measurements can complement 
each other and be of great usefulness. 

In nuclear astrophysics, the 8“decay half-lives set 
the timescale of the rapid neutron capture process (r- 
process), and hence influence the production of heavy 
elements beyond iron in the universe m3- Other weak 
processes like electron capture are of paramount impor¬ 
tance in the simulation of the latest stages of the evo¬ 
lution of massive stars and in the understanding of pos¬ 
sible supernova explosions. In particle physics, nuclear 
8-decay was used to obtain the first experimental ev¬ 
idence of parity-violation. The superallowed Fermi 8~ 
decay of nuclei can be exploited to verify the unitarity of 
the Cabibbo-Kobayashi-Maskawa matrix jsj, Q. 

Extensive studies of 8"decay have been carried out 
both experimentally and theoretically. In experiment, 
important advances in the measurement of the decay 
half-lives have been achieved in recent years with the 
development of radioactive ion-beam facilities. As for 
theory, many approaches have been formulated ranging 
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from gross theory El to the interacting boson-fermion 
model a Ab-initio approaches can still be used only 
for light nuclei 111 , whereas the nuclear shell-model can 
only cover the nuclear chart up to intermediate values of 
A « 40-50 and/or around magic regions if some frozen 
core is assumed. It performs quite well in reproducing 
8-decay half-lives in these cases (cf., e.g., the extensive 
studies in the sd shell [13], in the pf shell 0 , and in 
heavy nuclei 013 ; see also 0 and references therein). 

Another approach, which can be applied throughout 
the whole isotope chart, is based on the random phase 
approximation (RPA), and on its extension to superfluid 
systems, namely the quasi-particle RPA (QRPA). Many 
versions of QRPA have been applied for the study of 
8-decay, and most of them are based at least in part 
on phenomenological ingredients 0-(2l|. Discussing the 
pros and cons of the various phenomenological inputs is 
outside the scope of this work. Certainly, it would be 
desirable to reproduce the 8-decays half-lives within the 
framework of a self-consistent model without adjustable 
parameters. 

Actually, nuclear 8-decay has also been investigated 
within self-consistent (Q)RPA models based either on 
nonrelativistic or relativistic effective interactions (or en¬ 
ergy functionals); in these approaches the 8-decay half- 
lives are usually overestimated j22l - (25j . This deficiency 
can be to a good extent cured in the case of open- 
shell nuclei with the inclusion of an attractive isoscalar 
proton-neutron (pn) pairing force, that can shift some 
GT strength downwards and thus reduce the half-lives. 
However, the isoscalar pn pairing has no, or little, effect 
on closed- or subclosed-shell nuclei like for example 78 Ni 
and 132 Sn [2^ - l25| . Thus, the half-lives of such nuclei can¬ 
not be reduced; in some cases 132 Sn is even predicted to 
be stable El- Therefore, there must exist other corre¬ 
lations, other than isoscalar pn pairing, that are capa¬ 
ble to reduce such half-lives. One possible candidate is 
an attractive tensor force, as it has been pointed out in 
Ref. 0]. It must be stressed that the introduction of 
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new terms in the Hamiltonian, like isoscalar pn pairing 
or tensor terms, requires the tuning of one or more pa¬ 
rameters. It is of considerable interest to check whether 
the inclusion of new correlations without the fitting of 
new parameters can improve the agreement_of /3-decay 
half-lives with experiment. A recent study j27i| suggests 
that coupling of the GT states with two-phonon states, 
as well as the effect of the tensor force, increases the /3- 
decay rates. However, the model used in Ref. [27[ is not 
fully self-consistent since it adopts the Landau-Migdal 
approximation for the residual force, and makes a strong 
selection of two-phonon states. To investigate the role of 
correlations beyond mean-field in a self-consistent model 
and in a larger model space, we apply for the first time 
the RPA plus particle-vibration coupling (PVC) to the 
study of /3-decay. 

The RPA approach for giant resonances or other vibra¬ 
tional states is restricted to configurations of one particle- 
one hole (lp-lh) nature. In order to reproduce the ob¬ 
served spreading widths, one must consider the damping 
caused by the coupling to more complicated states, like 
2p-2h configurations [28)-[30|. An effective way to account 
for (most of) the observed spreading widths is to include 
the contribution of lp-lh -1 phonon configurations or, in 
other words, to take into account the coupling of single¬ 
nucleon states to the collective low-lying (mainly surface) 
nuclear vibrations or phonons [3lJ . We call this model 
RPA plus PVC, and we note that self-consistent versions 
of such a model have been realized based on relativis¬ 
tic and nonrelativistic interactions or energy function¬ 
als. Good agreement with the experiment has been ob¬ 
tained in the case of line shape of the GT and spin-dipole 
strength distribution l32M.'lol . In our application of the 
self-consistent RPA+PVC model, based on Skyrme en¬ 
ergy density functionals [34, HI, it was found that the 
coupling to phonons produces a downward shift of the 
GT excitation strength distribution (by about 1-2 MeV), 
accompanied by the development of broadening and/or 
fragmentation. In short, we expect that correlations be¬ 
yond lp-lh ones, and due to PVC, can play an important 
role also in the /3-decay processes. Accordingly, in this 
Letter, we apply our Skyrme RPA plus PVC model to 
the study of /3-decay. Our model is self-consistent in the 
sense that the same Skyrme force is used to calculate the 
single-particle levels and the RPA spectrum (both the 
GT one and those of the low-lying surface vibrations to 
be coupled to it), as well as the PVC vertices. 

The PVC effects are introduced along the line of Refs. 
[HI 1H1 as we shall now discuss. The RPA in its standard 
matrix form produces a set of 1 + eigenstates |n), having 
energies E n and strengths B n , as well as the forward¬ 
going and backward-going amplitudes denoted by X^ 

and , respectively. We then couple these RPA states 
with a set of doorway states consisting of a p-h excita¬ 
tion coupled to a collective vibration. These collective 
vibrations are obtained by computing the RPA response 
with the same Skyrme interaction. We have considered 
phonons with natural parity J n = 0 + , 1 _ , 2 + , 3 _ , 4 + , 
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FIG. 1: (Color online) /3-decay half-life of 78 Ni, calculated 
by RPA and RPA+PVC approaches with several different 
Skyrme interactions, in comparison with the experimental 
value [3§ |. 

5 _ , and 6 + having energies smaller than 20 MeV and as¬ 
sociated with a fraction of the total (isoscalar or isovec¬ 
tor) strength larger than 5%. 

The self-energy of the RPA state \n) is given by 

Y n {E M ) = ]T w£ hyh ,(E M )X$xW 

ph,p'h' 

+W^ h ,(-E M )Y^Y% f , ( 1 ) 

where the matrix elements w£ h p i h >(E M ) are spreading 
terms associated with the coupling of the lp-lh configu¬ 
rations with the doorway states, defined in Refs. [34j.l35|. 
They are complex and energy-dependent, calculated by 
using an averaging parameter p that avoids divergences 
and represents, in an approximate way, the coupling of 
the doorway states to even more complicated configu¬ 
rations. We have found that the lifetimes converge for 
small values of 77 , which are expected to be appropri¬ 
ate for low-lying, discrete states. In our calculations we 
have set r/ = 0.05 MeV. One can then calculate the GT 
strength distribution from the Gaussian averaging [36} 

1 (g a ,-g n -AE„) 2 

S{E M ) = £- 7 ^=e ^ B n , (2) 

„ 0+V27T 

where a n is defined as a n = (^ + 17 )/\/ 21 n 2 , with A E n = 
ReYj n (EM) and r„ = -2ImE n (Em)- 

Once the strength function has been obtained, the /3- 
decay half-life of an even-even nucleus is calculated in 
the allowed GT approximation by using the expression 

mil 


9a J AtiH S(E M )f(Z, Lo)dE M 

with D = 6163.4 s and gA = 1- The integration includes 
all final 1 + states having an excitation energy Em , re¬ 
ferred to the ground state of the mother nucleus, smaller 
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than A rlH = 0.782 MeV, which is the mass difference 
between neutron and hydrogen. If the energy is instead 
referred to the ground state of the daughter nucleus, and 
is denoted by E, one has E = Em — AB, AB being the 
experimental binding energy difference Bm — Bp. This 
choice is often convenient, because the calculated energy 
of the final 1 + states can then be directly compared to 
the experimental spectrum of the final nucleus. It may 
happen that the calculated energy of the lowest state lies 
at negative energy, that is, below the experimental value 
of the ground state energy of the daughter nucleus. The 
upper limit of integration in (£JJ becomes equal to the Qp 
value (A n H — A B = Qp), namely to the atomic mass dif¬ 
ference between the mother and daughter nucleus. Eq. 
d3]) then becomes 


9a f Qp S(E)f(Z, uj)dE' 

The integrated phase volume f(Z,ui) is given by 

f‘UJ 

f(Z,u)= / p e E e (co-E e ) 2 F 0 (Z + l,E e )dE e , (5) 

J m e c 2 

where p e ,E e , and Fq(Z + 1 ,E e ) denote the momentum, 
energy and Fermi function of the emitted electron, re¬ 
spectively. u) is the energy difference between the ini¬ 
tial and final nuclear state, connected with GT energy 
E (or Em) by ui = Qp + m e c 2 — E = A np — Em where 
A„ p = 1.293 MeV is the mass difference between neutron 
and proton. 

Fig. |T| shows the results for the /3-decay half-life of 
78 Ni, a doubly magic nucleus that is an important waiting 
point in the r-process, and represents as a major bottle¬ 
neck in the synthesis of heavier elements. As recalled in 
the beginning of this Letter and evident in the figure, the 
different Skyrme interactions are not well constrained in 
the spin-isospin channel and their predictions for the half- 
life at the RPA level can vary over more than two orders 
of magnitude. For instance, SAMi f3§j predicts a rather 
collective GTR and does not leave much strength at low 
energy, resulting in a long half-life, at variance with SLy5 
g3| that has weak spin-isospin terms and produce a less 
collective GTR, leaving more strength at low energy and 
providing a much shorter half-life. The parameter sets 
SAMi, SGII jU, SkM* [42j and Sill j+| overestimate 
the half-life, while the interactions SLy5 and Skx jUJ are 
in agreement with the experiment at the RPA level. 

The inclusion of PVC effects reduces the half-lives for 
all interactions systematically. The reduction factor R is 
larger for SAMi ( R sa 42) and SGII (R ss 10) while it is 
equal to about 4 for the other four interactions. Within 
RPA+PVC, the half-life obtained with the sets SkM* 
and Sill falls within the experimental error. It has to 
be stressed that the Skyrme force SkM* does not only 
reproduce well this /3-decay half-life, but also describes 
well the giant resonance line shape in 208 Pb and 56 Ni at 
the PVC level (H), IH]- We will show results for the /3- 
decay of 132 Sn, 68,78 Ni, and 34 Si with all interactions, but 
we will discuss in detail the SkM* results. 


In order to understand the reasons for the systematic 
decrease of the half-lives after the inclusion of phonon 
coupling, we display in Fig. [2] the GT strength dis¬ 
tributions (with respect to the daughter nucleus), the 
cumulative sums of the strengths, and the cumulative 
sums of 1/7 i/2 [that is, the values obtained from Eq. Q 
when Qp is replaced by a running E in the upper limit 
of the integral in the denominator]. Generally speak¬ 
ing, for all nuclei we have considered the GT peaks are 
shifted downwards when going from RPA to RPA+PVC. 
As seen in Fig. HT a), in the case of 132 Sn, in RPA there 
is no strength below Qp, so that the nucleus is stable. 
The lowest RPA state located at E = 3.6 MeV is moved 
into the /3-decay window and becomes very close to one 
of the experimentally observed 1 + states (which, how¬ 
ever, has a decay branching ratio of only I = 0.87%). 
The total strength below 3.6 MeV then does not change 
much from RPA to RPA+PVC [Fig. Elb)]. Another 
(lower) 1 + state, which experimentally has by far the 
largest decay branching ratio, is not reproduced by our 
model; thus, our RPA+PVC can predict a finite value 
of the half-life, providing a qualitative improvement with 
respect to RPA but still overestimates the experimental 
finding [Fig. Ufc)]. In the case of 68 Ni, RPA predicts a 
state in the /3-decay window but its energy is higher than 
experiment as seen in Fig. Hid). This explains why the 
half-life is overestimated with respect to the experimen¬ 
tal finding: the contribution of this state to I/T 1/2 is so 
small that it is only visible in Fig. Elf) when multiplied 
by a factor 10. With the inclusion of PVC, the RPA peak 
at 1.5 MeV is moved even slightly below the experimen¬ 
tal ground state energy. Although, as in the previous 
case, the strength of this state is not changed much as 
compared with RPA [Fig. H](e)], it gives a very large 
contribution to I/T 1/2 because of the large phase-space 
factor. As a consequence, the half-life is smaller than in 
experiment. In the case of 34 Si, in RPA one finds three 
peaks located at E = —0.86, 3.1, and 4.2 MeV. The first 
one lies below the experimental ground state and deter¬ 
mines the value of 1 /T 1/2 [Fig- Eli)]- This peak carries a 
very small value of the strength and therefore the exper¬ 
imental lifetime is largely overestimated. Including the 
PVC the strength becomes fragmented [Fig. Elg)]- One 
can identify five peaks at E = —2.2,1.0,1.7, 2.6, and 3.1 
MeV, contributing respectively 15%, 49%, 24%, 3%, and 
9% of the total value of I/T 1 / 2 , which becomes much 
larger than in RPA, substantially improving the agree¬ 
ment with the experimental lifetime. For the nucleus 
,8 Ni, the small strength at E = 5.6 MeV gives almost 
all the contribution to I/T 1/2 in the RPA model [Fig. 
EH)] which underestimates the experimental value. With 
PVC, the state at E = 5.6 MeV is shifted to a state 
peaked at 4.0 MeV [Fig. EI.j)] with strength unchanged 
[Fig. Elk)]. As a result, due to the increase of phase 
space its contribution to I/T 1/2 becomes about 3.4 times 
larger than in RPA [Fig. Ell)]- The strength distribution 
above this peak contributes 22% of the total l/Ti/ 2 - 

The resulting calculated lifetimes for these four nuclei, 



4 


_ 40 
\ 30 
~ 20 
<5 10 
« o 

_ 3 
O 2 

w 1 
0 

« 0.2 

o 

" 0.1 
t 

"■ 0.0 



oo 

a 


FIG. 2: (Color online) Experimental data related to /?-decay from nuclei 132 Sn, 68 Ni, 34 Si, and 78 Ni are compared with 
theoretical results obtained with the SkM* interaction. In these panels, the excitation energies Em calculated with respect 
to the mother nucleus are transformed to E, the excitation energies referred to the ground state of daughter nucleus, using 
experimental binding energy difference (see the text); accordingly, the vertical dotted lines show the experimental value of 
Q t 3 [45|. Top panels: GT _ low-lying strength associated with the discrete RPA peaks _B(GT“) (dashed lines) and with the 
continuous RPA+PVC strength distributions 5(GT“) (solid lines). The arrows indicate the experimental energies of the 
measured 1 + states [45]. Middle panels: cumulative sum of the RPA and RPA+PVC strength shown in the top panels. Bottom 
panels: cumulative sum of l/Ti/ 2 - The experimental values of 1 /Ti/ 2 01] for each nucleus are indicated by the stars. The 
strength of the lowest RPA and RPA+PVC peaks in panel (g) and the RPA l/Ti /2 in panel (f) have been multiplied by a 
factor of 100 and 10, respectively. 



FIG. 3: (Color online) The /3-decay half-lives of 132 Sn, 
68 Ni, 34 Si, and 78 Ni, calculated by RPA and RPA+PVC ap¬ 
proaches, respectively, in comparison with experimental val¬ 
ues [Iff]. The arrows denote half-lives longer than 10 6 s. 


in the case of all the forces used in Fig. [T) are compared 


with experiment in Fig. [3] The RPA results generally 
markedly overestimate the half-lifes for all nuclei. An 
exception is represented by the interaction Skx, in which 
case one obtains a good agreement with data at the RPA 
level; this is associated with the fact that the properties of 
132 Sn, 68 Ni, and 34 Si, as well as the single-particle levels 
of 132 Sn and 34 Si, have been used to fit of the parame¬ 
ters of this force [H[. The effect of the PVC decreases 
the values of TV 2 by large factors compared to RPA, 
substantially improving the agreement with experimen¬ 
tal data, except for Skx and (partially) for SLy5. With 
the inclusion of PVC effect, the interactions SkM* and 
Sill give the best agreement with data. More in detail, 
in the case of SkM*, the lifetime is still large in 132 Sn and 
small in 68 Ni, in keeping with the errors in the position 
of the lowest 1 + state (cf. Fig. [2]). Theory agrees instead 
very well with data in the case of 34 Si and ' 8 Ni. 

In conclusion, we have shown that, starting from RPA, 
the coupling between particles and vibrations causes a 
significant downward shift in the GT strength function 
of the four nuclei 132 Sn, 68 Ni, 34 Si, and 78 Ni (treated as 
magic). The /3-decay half-life is more sensitive to the po¬ 
sition of the 1 + states rather than to the strength, which 
is not much changed in going from RPA to RPA+PVC. 
This is due to the strong increase of the decay phase 
space factor as the energy decreases. As a consequence, 
the lifetime is reduced in the case of RPA+PVC and 
the agreement between theory and experiment is in gen- 


















































5 


eral substantially improved. In particular, the interaction 
SkM* that had been previously shown to perform well in 
magic nuclei as far as the line shape of the GT resonance 
is concerned [35|, leads to overall good agreement with 
/1-decay data. 

We can expect that including the effect of PVC will 
also be helpful in the case of other weak interaction pro¬ 
cesses. For example, in the case of the electron capture 
on fp-shell nuclei, the threshold position and cross sec¬ 
tions at low energies are underestimated with respect to 
e.g. shell-model calculations [I(|. PVC can lower the 


excitation energies in the GT + channel in a similar way 
as it has been shown here. The study of open-shell nu¬ 
clei by including pairing correlations is envisaged. Then 
the model can be employed to predict the half-lives of 
r-process bottle neck nuclei with N = 82, which play 
an important role for the duration of the r-process, and 
hence can help to understand the origin of heavy elements 
in the universe. 
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